Direct observation of the superconducting gap in phonon spectra 
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We show that the superconducting energy gap A can be directly observed in phonon spectra, as 
predicted by recent theories. In addition, since each phonon probes the gap on only a small part 
of the Fermi surface, the gap anisotropy can be studied in detail. Our neutron scattering investiga- 
tion of the anisotropic conventional superconductor YNi2B2C demonstrates this new application of 
phonon spectroscopy. 

PACS numbers: 74.70.-b, 74.72.-h, 63.20.dd, 63.20.kd, 78.70.Nx 



In their seminal paper [lj, Axe and Shirane found that 
the onset of sueprconductivity in NbaSn affects the fre- 
quencies and the linewidths of phonons whose energies 
are comparable to the superconducting energy gap 2A. 
Similar effects were observed by neutron and Raman scat- 
tering in a number of conventional superconductors 0, 
but also in high-T c copper oxides [1, 0, H, H, ■ in many 
of these cases the authors argued that it was possible to 
determine whether the phonon energy was above or be- 
low 2A from the comparison of phonon frequences and 
linewidths above and below T c . Here we show that in 
the case of strong electron-phonon coupling, the entire 
phonon lineshape can change dramatically accross T c al- 
lowing a precise determination of 2 A. 

A few years ago, strong superconductivity-induced 
changes of the phonon lineshape were reported for a par- 
ticular phonon in YNi2B2C[8j and LuNi2B2C[9], which 
- for the first time - could not be described by a simple 
change in the phonon linewidth. This observation stim- 
ulated further theoretical workfiol [Tl| . We performed 
very precise inelastic neutron scattering measurements of 
many different phonons in YNI2B2C (T c = 15.2 K) with 
a much smaller statistical error covering a much larger 
region of reciprocal space than in Ref. [8|. 

The neutron scattering experiments were performed on 
the IT triple-axis spectrometer at the ORPHEE reactor 
at LLB, Saclay, and on the PUMA triple-axis spectrom- 
eter at the research reactor FRM II in Munich. Double 
focusing pyrolytic graphite monochromators and analy- 
sers were employed in both cases. A fixed analyser en- 
ergy of 14.7 meV allowed us to use a graphite filter in the 
scattered beam to suppress higher orders. The experi- 
mental resolution was obtained from calculations using 
a formalism checked on phonons with very low coupling 
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strength. The wave vectors are given in reciprocal lat- 
tice units of (2n/a 2ir/b 2-k/c), where a = b = 3.51 A 
and c = 10.53 A. The single crystal sample of YNi2B2C 
weighing 2.26 g was grown by the floating zone method 
using the n B isotope to avoid strong neutron absorption. 
The single crystal sample was mounted in a standard or- 
ange cryostat at LLB and in a closed-cycle refrigerator 
at FRM II, allowing us measurements down to T = 1.6 K 
and 3K, respectively. 

Figure [1] shows the evolution through T c of the low 
temperature lineshape of the phonon at the endpoint 
of the transverse acoustic branch in the (110) direction. 
This phonon has never been investigated before but its 
strong coupling to electrons has been predicted by den- 
sity functional theory [l^. It appears at q = (0.5, 0.5, 
0) (the M-point). On cooling from T = 200 K to a tem- 
perature just above T c — 15.2 K, this phonon softens by 
about 10 % and broadens substantially, indicative of a 
strong electron-phonon coupling (data not shown). How- 
ever, the lineshape remains Lorentzian to a very good 
approximation [18j. On further cooling through T c , the 
lineshape starts to deviate strongly from a Lorentzian. 
In particular, a step-like increase at a certain energy E s 
appears as the temperature is lowered through T c (fig. [TJ 
top). E s increases with decreasing temperature below 
T n . According to the theory developed by Allen et al. 
[lfj |. a part of the low energy tail which lies below the 
gap, is pushed up in energy to form a narrow spike at 
2 A (fig. [1] bottom). The theory is based on the full quan- 
tum mechanical treatment of electron-phonon coupling 
where vibrational and electronic excitations mix into hy- 
brid modes. The finite spectrometer resolution should 
wash out the spike in experiments resulting in a sharp in- 
tensity increase at 2A. The predicted intensity increase 
corresponds to the one we observe at E s . We show below 
that other experimental features are also in an excellent 
agreement with this scenario. Thus we will interpret our 
results using the theoretical framework developed in Ref. 
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A close comparison of the calculated lineshape with the 
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FIG. 1: (color online) top: Evolution of the neutron- 
scattering profile measured on YM2B2C at Q = (0.5,0.5,7) 
at T = T C (15.2K) and below T c . Solid lines are guides to 
the eye. The broken line denotes the background as deter- 
mined from measurements at neighboring q-points. bottom: 
Calculated phonon lineshapes based on the theory of Allen et 
al.fiol]. using parameters extracted from the experimentally 
observed lineshape in the normal state (thin lines). Thick 
lines are obtained after convolution with the experimental 
energy resolution. 



low energy part of the observed phonon spectral function 
can give highly accurate values of 2 A as a function of tem- 
perature (fig. [2 left). Close to T c , where the step height 
is small, it is more convenient to scan the temperature at 
a fixed energy transfer: a sudden intensity decrease on 
cooling will mark the temperature where the gap is equal 
to the energy transfer (fig. [2 right). 

The theory of Allen et al.[l(| predicts not only the 
suppression of the spectral weight below 2A just men- 
tioned, but also a downward shift of the intensity max- 
imum and an intensity build-up on the low energy side 
of the peak (fig. [fl bottom) . The observed changes are 
qualitatively similar, but are somewhat stronger than 
predicted. These discrepancies cannot be remedied by 
tuning the three parameters of the calculations. Two of 
them, the phonon energy and the phonon linewidth ob- 
served just above T c , are entirely fixed by the phonon 
properties in the normal state and the third, the gap 
value is practically fixed as well by the sharp intensity 
increase at E s . Notwithstanding the shortcomings of the 
actual theory - which might be due to the weak coupling 
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FIG. 2: (color online) left: Background subtracted energy 
scans taken at Q = (0.5,0.5,7) at three different temperatures 
including an offset of 300 counts for clarity. Lines were cal- 
culated from the theory of Allen et al.[lfj|] after scaling to the 
experimentally observed step-heights and folding with the ex- 
perimental energy resolution. The vertical bars indicate the 
gap values used in the calculations, respectively, right: Tem- 
perature scans taken at Q = (0.5,0.5,7) for three different 
energy transfers including an offset of 35 counts for clarity. 
Lines are as in the left panel. Vertical bars indicate the tem- 
perature for which the energy transfer is equal to the gap 
value 2A in the calculation. 



approximation -, it predicts the observed changes with 
temperature over a wide range of energies remarkably 
well (fig.©. 

The features in the intensity vs. temperature curves 
(fig.[H right) related to the gradual opening of the gap are 
quite sharp and, therefore, the gap values extracted from 
the data depend very little on fine details of the model. 
The temperature dependent gap value deduced in this 
way is plotted in fig. 0J The experimental curve deviates 
somewhat from the BCS weak coupling curve, which is 
not surprising for a superconductor with a T c as high as 
15.2 K. Likewise, the value of the zero temperature gap 
indicates strong coupling because the ratio 2A//csT c = 
4.3 (kg is the Boltzmann constant) exceeds considerably 
the weak coupling limit of 3.52. 

The method of extracting the superconducting gap 
from phonon lineshapes explained above looks straight- 
forward. Therefore, one might ask why it has not yet 
been used so far. There are two stringent conditions for 
making the replica of the superconducting gap clearly 
visible: the phonon energy has to be considerably larger 
than 2A but on the other hand, the linewidth in the 
normal state has to be large enough that the tail of the 
phonon line extending into the gap region can still be dis- 
criminated from the background. One might also think 
to study phonons with an energy comparable to the gap 
energy. Here, the superconductivity-induced changes of 
the phonon lineshape are fairly strong, but extraction of 
the gap energy is possible only with the help of detailed 
calculations. An example for such a case is the acoustic 
phonon in the (lOO)-direction at wave vector q = (0.5, 
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FIG. 4: (color online) Superconducting energy gap 2A versus 
temperature as evaluated from phonon data. Broken lines 
are calculated from weak-coupling BCS theory scaled to the 
experimental values at T = 3K. Red and blue points were 
deduced from phonon measurements at different momentum 
transfers. Lines are guides to the eye. 



FIG. 3: (color online) top: Evolution of the neutron- 
scattering profile measured on YN12B2C at Q = (0.5,0,8) 
above and below T c . Solid lines are guides to the eye. The bro- 
ken line denotes the background as determined from measure- 
ments at neighboring q-points. bottom: Calculated phonon 
lineshapes based on the theory of Allen et al.[l(il]. using pa- 
rameters extracted from the lineshape observed in the normal 
state (thin lines). The thick lines are obtained after convolu- 
tion with the experimental resolution. 



0, 0). Profound changes of the lineshape of this phonon 
were first reported by Kawano et al.[8[. Puzzled by the 
appearance of an intense and relatively narrow line be- 
low T c Kawano et al. thought to have found a new sharp 
excitation. We reinvestigated this phonon and show that 
the peculiar lineshapes appearing below T c can be quite 
naturally understood using again the theory of Allen et 
al.pil. 

The phonon at q = (0.5, 0, 0) has a significantly lower 
energy than the M-point phonon discussed above (7meV 
as compared to 10 meV) and moreover, has a much larger 
linewidthpjl. As a consequence, the superconductivity- 
induced redistribution of spectral weight is even stronger, 
but the superconducting gap cannot be inferred from the 
data as easily, except for temperatures close to T c where 
the gap is still small. For lower temperatures, most of 
the spectral weight condenses into a fairly sharp peak 
whose energy is somewhat below the gap. According to 
Ref. Qil this resonance can be regarded as a mixed vibra- 
tional/superelectronic collective excitation. Its position 
with respect to 2A has to be determined from calcula- 
tions for the parameter values of this particular phonon. 
Again, the theory does not reproduce the observed line- 



shapes quantitatively, but sufficiently well to precisely 
determine the gap value (fig. [3]). 

Previous publications [HQ observed line narrowing for 
phonons with an energy smaller than 2A. We found that 
this temperature dependence can be quantitatively ac- 
counted for by the theory of Allen et al.fioj. Unfortu- 
nately, there are no phonons in YNI2B2C with energies 
< 2 A, which have a sufficiently strong coupling to search 
for temperature effects of the linewidth in this compound. 

The low temperature lineshapes of the two phonons 
discussed so far look quite distinct. A study of several 
phonons along the (lOO)-direction reveals that the line- 
shape may take on any intermediate form, depending on 
the phonon energy relative to the gap and on the normal 
state phonon linewidth. When going from q = (0.5, 0, 0) 
towards the zone boundary, the normal state phonon en- 
ergy increases and the linewidth shrinks. As a result, the 
low temperature lineshape gradually evolves from that 
observed at q = (0.5, 0, 0) to that of the M-point phonon 
(fig. OH (a)-(c)). Fig. [5] is interesting also for another 
reason: the superconductivity-induced redistribution of 
spectral weight involves the optic line at 14 meV as well 
with the spectral weight near both peaks transfered to 
much lower energies near the gap. At these q-points, the 
acoustic and the optic modes strongly hybridize. Un- 
derstanding this phenomenon in detail would require a 
highly nontrivial extension of the theory, which has not 
yet been made. 

When plotting 2A extracted from the various phonons 
studied by us, it becomes clear that they probe different 
energy gaps, ranging from 4 meV to more than 6 meV 
(figs. 4, 5(d)). Phonons probe the gap on the parts of 
the Fermi surface connected by the phonon wave vector. 



4 



100 







(a) (0.65,0,8) 

I \f X 


(b) (O./'.O.SS 

• 17 K O 

o 2 K F% 






(c) (0.75,0,8) 

* J, S 


4 


8 12 16 
E (meV) 


_ 6 
> 

1 < 

CM 

4 


* 

<f * T = 2K 





4 8 12 16 0.4 0.6 0.8 

E (meV) (x, 0, 0) (r.l.u.) 



FIG. 5: (color online) Panels (a)-(c) show background sub- 
tracted energy scans taken at Q = (0.65-0.75,0,8) at T = 2K 
(open symbols) and 17 K (filled symbols). Lines are guides to 
the eye. Panel (d) describes the evolution of the low tempera- 
ture superconducting gap value extracted from phonon scans 
along the (100) direction. 



Only theory can tell which parts are actually involved. 
Density functional theory calculations may provide such 
information, but they require the technically involved so- 
lution of the anisotropic gap equations including the full 
momentum dependence of the electron-phonon coupling, 
which has not yet been carried out for YNi2B2C. In spe- 
cial cases, however, the respective parts of the Fermi sur- 
face can be inferred from the Fermi surface topology: for 
instance, the extremely large linewidth of the q = (0.5, 
0, 0) phonon can be attributed to a nesting phenomenon 
predicted by theory [T^- 



Whereas the intensity step observed at q = (0.5, 0.5, 
0) is very sharp, i.e. resolution- limited (» lmeV), it is 
significantly broader (w 2mcV) at q = (0.7, 0, 0) and 
(0.75, 0, 0), and therefore cannot be linked to a single, 
well-defined gap value. In this case the phonon wave 
vector may connect several different pieces of the Fermi 
surface probing a superposition of gap values. 

Our results for the gap anisotropy in YM2B2C fit 
nicely to those extracted from tunnelling data using 
point-contact spectroscopy [l3j. Further, the lowest gap 
seen in our phonon experiments agrees very well with 
the value derived from critical field H c2 (T) data[3|- 
On the other hand, we have no evidence that the gap 
goes below 2 A = 4meV anywhere on the Fermi surface 
and possibly even to zero, as was proposed on the basis 
of field-angle-dependent heat capacity measurement sllal 
and field-angle dependent thermal conductivity data[16(. 
Extensive energy and wave vector resolved density func- 
tional theory calculations I12f , which are in very good 
agreement with experiment (1 7| . did not show any other 
suitable phonons in order to further probe the possible 
presence of nodes. 

To conclude, we confirm the prediction of Allen et al. 
that the phonon line may acquire very peculiar shapes 
with prominent features directly related to the super- 
conducting energy gap 2A. In contrast to a theoretical 
proposition [llj, this effect is not restricted to phonons 
with wave vectors close to an extremum vector of the 
Fermi surface but may appear throughout the Brillouin 
zone. The gap value 2A can be directly inferred from the 
lincshape, often without recourse to a fit. This allows ex- 
traction of the gap values as a function of temperature 
with high accuracy. By studying phonons of different 
wave vectors, even the gap anisotropy can be explored 
by phonon spectroscopy. Thus phonon spectroscopy can 
be a much richer source of information on the supercon- 
ducting properties than hitherto assumed. 
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